arXiv:1509.00637vl [math.PR] 2 Sep 2015 


Clustering of high values in random fields 


Helena Ferreira* Lui'sa Pereiraf Ana Paula Martins^ 


Department of Mathematics, University of Beira Interior, Portugal 

September 3, 2015 


Abstract: The asymptotic results that underlie applications of extreme random fields often 
assume that the variables are located on a regular discrete grid, identified with 7 ?, and that 
they satisfy stationarity and isotropy conditions. Here we extend the existing theory, con¬ 
cerning the asymptotic behavior of the maximum and the extremal index, to non-stationary 
and anisotropic random fields, defined over discrete subsets of We show that, under a 
suitable coordinatewise long range dependence condition, the maximum may be regarded as 
the maximum of an approximately independent sequence of submaxima, although there may 
be high local dependence leading to clustering of high values. Under restrictions on the lo¬ 
cal path behavior of high values, criteria are given for the existence and value of the spatial 
extremal index which plays a key role in determining the cluster sizes and quantifying the 
strength of dependence between exceedances of high levels.The general theory is applied to 
the class of max-stable random fields, for which the extremal index is obtained as a function 
of well-known tail dependence measures found in the literature, leading to a simple estima¬ 
tion method for this parameter. The results are illustrated with non-stationary Gaussian and 
1-dependent random fields. For the latter, a simulation and estimation study is performed. 
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1 Introduction 

Extremes of variables like wind, temperature and precipitation can affect anybody at any place. The 
potential consequences include increases in severe windstorms, flooding, wildfires, crop failure, population 
displacements and increased mortality. Apart from their direct impacts, these events will also have 
indirect effects such as increased costs for strengthening infrastructure or higher insurance premiums. 
When the interest lies in the study of variables measured at specifically-located monitors, such as the 
variables mentioned above, as well as air pollution, soil porosity or hydraulic conductivity, among others, 
spatial modeling is necessary, so random fields constitute an active area of current research. 

The treatment of spatial and temporal dependence in random fields has been influenced by the mul¬ 
tivariate Gaussian model, where the dependence is characterized by the covariance structures. However, 
this model excludes all the situations of marginal distributions with heavier tails than the Gaussian dis¬ 
tribution, leaving aside a huge set of problems related to rare events. Extreme Value Theory plays an 
important role in these situations. 
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A considerable amount of work has been done in extending results of Extreme Value Theory to 
random fields which have 7? as their parameter space. Although their lack of easy separation of past 
and future, a general version of the classical Extreme Types Theorem was given and the existence of the 
extremal index shown, by replacing a single global dependence restriction by several assumptions, each 
dealing with one coordinate direction, for which past-future separation is considered m. m , among 
others). Under local restrictions on the oscillations of the values of the random field, Ferreira and Pereira 
([S]) and Pereira and Ferreira (|I31) compute the extremal index from the joint distribution of a finite 
number of variables. 

In a random field with high local dependence, an exceedance is likely to have neighboring exceedances, 
resulting in a clustering of exceedances, which leads to a compounding of events in the limiting point 
process of exceedances m)- 

The aforementioned results assumed that the variables are located on a regular grid, identified with 
Z^, and sometimes that they satisfy stationarity and isotropy conditions. This is a big restriction for the 
majority of the applications since usually spatial data are not regularly spaced, stationary and dependence 
is anisotropic, due to the presence of a main direction of dependence. 

In this paper we extend the existing theory, concerning the asymptotic behavior of the maximum, 
to non-stationary and anisotropic random fields, Z 5 = {Z{x) : x G S'}, where S = Un>i 
^ = {^n}n>l is an increasing sequence of sets of isolated points of subject to conditions on long 
range and local dependencies. We will assume, without loss of generality, that the variables Z{x), x £ S, 
have common distribution F, being F the corresponding survival function. We will denote the maximum 
and the minimum of Z{x) over B C S hy \/^^^Z{x) and /\^^qZ{x), respectively. More precisely, in 
Section 2 we define an asymptotically independence condition under which we prove that 
n > 1 , behaves asymptotically as independent maxima over a family of disjoint subsets of An- 

The way spatial extreme events interact is also of interest in spatial statistics. For example, an 
unusually stormy day at a particular location may be followed by another one at the same or a neighboring 
location. This type of dependence among spatial extremes can be summarized through the spatial 
extremal index of the sequence = {Z{x) : x £ A„}„>i. 

Definition 1.1. The sequence Za has spatial extremal index 0 a if, for each r > 0 and any sequence of 
real numbers {'iin('r)}„>i satisfying 


\xG.A7t, / 

where Ua denotes the indicator function of the event A, it holds that 


( 1 . 1 ) 


lim pi \/ Z(x) <'Un(r) ) = exp(-0Ar). 

t^ + OO \ '' / 


The extremal index of Za is the key parameter to relate the limiting distributions of \/x&An 
and VisAn E{x), where Za = {Z{x),x £ An}n>i is a sequence of independent and identically random 
variables having the same distribution function F as each variable of the sequence Za- In fact, if f{n) is 
the number of locations on An and there exists a sequence of real numbers {un{T)}n>i satisfying ( 1 . 1 ), 
then 


P 


\J Z(x) < Unix) j = FE'^){un{T)) 

CG An / 


n^+oo 


exp 


- lim f{n)F{un{T)) = e~ 

n^+oo 


and 


P ( \l Z(x) < Unix) 

\xeA„ / 


> + 00 


exp 


- lim 9f{n)F{un{x)) = e 

n^+oo 


— ^-Sat 
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So, 


1. lim p[ \l Z{x) < ttn(r) ) = lim P [ \/ Z{x) < Unir) ) , with tJPn = 0Af{n) < '^An, that 

n—^+oo \ ’ / n—>+oo \ ' / 

\x&An / \xGBn / 

is, for the sequence of real levels {un{T)}n>i, \/ Z(x) behaves asymptotically as the maximum 

X&An 

of less than f{n) independent variables. 


2. ^liin^P( Y Z{x) < Unix)) = ^liin^P( \J Z{x) < Vn{0r) where VniOr) = F ^ > 

xeAn xeA„ 

Unix) = F~^ ^1 — that is, \J Z{X) behaves asymptotically as the maximum of the same 

x^A^ 

number of independent variables but relatively to a level higher than Un(x). 

We may then deduce that the limit of the sequence {cn = P (Va;eA„ — Un(x))}n>i is the same 
as the one we would obtain when considering {cn ^ P if iii we replace 

the levels Un{x), n > 1, by ’’appropriately close” levels Vn{x'), n > 1, with x' < x, or if we consider a 
sequence {Pn}n>i with Bn C An and jjPn ~ 0f{n), instead of sequence {An}n>i- This suggests that 
over An and considering the levels Un{x) we should not obtain isolated exceedances of Z{x) contrarily to 
Z{x), and therefore in this situation they occur in clusters. Later on we will prove a result that reinforces 
this intuition. 

We finish Section 2 with an existence criteria for the extremal index of Z^. 

Section 3 contains the theory surrounding the maximum and the extremal index of under re¬ 
strictions on its exceedance local path behavior, which allow clustering of high values. Surprisingly we 
obtain a simple method for computing the extremal index of sequence Z^ as the limit of a sequence of 
tail dependence coefficients. 

Section 4 is devoted to the application of the results to max-stable random fields. There we introduce 
the notions of local and regional extremal indices and relate them with 8a- Based upon these relations a 
simple estimator for 6a is given and its performance is analyzed with an anisotropic and non-stationary 
1-dependent max-stable random field. Conclusions are drawn in Section 5 and the proofs are collected 
in the appendices. 


2 Asymptotic spatial independence 

In this section, we show that, under a suitable long range dependence condition, the maximum of 
random fields defined over discrete subsets of M^, may be regarded as the maximum of an approxi¬ 
mately independent sequence of submaxima, even though there may be high local dependence leading to 
clustering of high values. 

The results are obtained through an extension of the methodology in Ferreira and Pereira l|13jL for 
extremes on a regular grid, relying on the novelty of irregularly occurring extremes in space. 

The dependence structure used here is a coordinatewise long range dependence condition, which 
restricts dependence by limiting | F (VxeCuD -^(^) - ^ (VxeC ■^(®) - P (Vxer> ■^(®) - | 

with the two index sets C,D C An being ’’separated” from each other by a certain distance In along each 
direction. 

Throughout we shall say that the pair (/, J) is in S{'Ki{An)-,ln) if / C 5 and J C S are subsets of 
consecutive values of Xi^An) separated by at least In values of Xi{An), where Xi, i = 1,2, denote the 
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cartesian projections. The cardinality of the sets and 7ri(An), ^ = 1,2, will be denoted by = f{n)-, 
'i,TTi{An) = fi{n), i = 1 , 2 and we will assume that /(n) —)• +oo as n —)■ +oo. 


Definition 2.1. Let {un}n>i be a sequence of real numbers. If there exist sequences of positive integers 
I = {ln}n>l ^ = {kn}n>l SUch that 


fi{n) 


In “hoc, kn ^ “hOO, knln r/ \ ' 

n^+oo n^+oo J yTl) n^+oo 


0 , for each i = 1,2, 


( 2 . 2 ) 


and kfa {In, Un) n^oo 

a{ln,Un) = sup 


pf V Z{x)<Un] -ply Z{x)<Un]pl\/ Z{x) 
\x&CUD / Vxec / \x&D 


< Ur, 


where the supremum is taken over sets C and D such that, for each i = 1,2, C, D £ S{'Ki{An),ln), we 
say that the sequence Za satisfies condition D{un,kn,ln)- 


Under D{un, kn, /n)-condition we have asymptotic independence of maxima over disjoint sets of loca¬ 
tions, as shown in the following result. 


Lemma 2.1. Suppose that the sequence Za satisfies condition D{un, kn, In) for a sequence of real numbers 
{Un}n>l such that 

{f{n)F{un)}n>i is bounded. (2.3) 

If I s,t = l,...,kn, are disjoint subsets of 7ri{An) and jW, s,t = l,...,kn, are disjoint subsets of 


T^ 2 {An), and 


= T^2 nVTi ^(/(®)) <r^An, s,t = l, ...,kn, 


we 

have 




( \ 


( \ 

P 

y Z{x) < Un 

-n^l 

y Z{x) < Un 1 


\xe / 

s,t 



0 . 


The next result proves that asymptotically the distribution of the maximum of Z over An, n > 1, 
coincides with the distribution of the maximum of Z over a union of conveniently chosen disjoint subsets 
of An, whenever condition D{un,kn,ln) holds. 

The underlying idea to obtain the asymptotic distribution of the maximum of Z over An, n > 1 > 1, 
is to subdivide An into kn disjoint subsets, Bn’^\ s,t = 1, ..., kn, using the following construction method 
of the family : s,t = 1 ,..., fenj: 

• build 7rj“^(/(^^) n An, •5 = 1, ...kn, with I^^\ 5 = 1, •••, kn, abutting subsets of consecutive values of 
7ri(yl„), maximaly chosen for the condition 

P{Z{x) > Un) < ^ ^ P{Z{x) > Un) 

XGtT^ (7^'^)) Pi An X^An 


• for each s = l,.../c„, build 7r2 n vr^ ^(/(^)) n An, t = l,...,kn, with t = l,...,kn, 

contiguous subsets of 7 r 2 ( 7 rj“^(/*^®)) n An) and maximally chosen such that 

P(Z(x)>U„)<^^P(Z(x)>Un). 

a;g7r^^( fj 7rj”^ (/(“)) p A„ ^ x£An 
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HiW 


Figure 1: Example of a set of disjoint blocks Bn’^\ s,t = 1, kn, with |J^ ^ Bn’^^ = An, for a particular 
An, n €N. 


Figure 1 llustrates one possible set of disjoint blocks Bn’^\ s,t = l,...,kn, with ~ 

constructed through the previous method, for a particular An, n G N. 


Lemma 2.2. If the sequence satisfies condition D{un,kn,ln), with {un}n>i verifying (2.3) then, for 
each n, there exists a family Bn of k") disjoint subsets of An, Bn'^\ s,t = l,...,kn, with ~ 

and such that 

( 


P ( Y Z{x) < tin I 

\x£A„ ) 


-P 


y Z(x) < Ur, 


0 . 


n—>-+oo 


totic independence of maxima over Bn'^\ s,t = 1 , 2 ,..., kn- 

Proposition 2.1. If the sequence satisfies condition D{un,kn,ln) with {un}n>i verifying (2.3), then, 
for each n, there exists a family Bn of kn disjoint subsets of An, Bn'^\ s,t = 1,..., kn, with ~ 

and such that 


As a consequence of Lemmas 2.1 and 2.2 we can now state the following result concerning the asymp- 


p( V Z{x)<Un]-llpl V 

\x&An / S,t 


Z{x) < Ur 


0 . 


n^+oo 


The following result gives a convenient existence criteria for the extremal index of Za and follows 

(s t) 

from Proposition 2.1: it depends on the local behavior of exceedances over Bn ’a s,t = 1,..., kn, namely, 
on the limiting mean number of exceedances of Un by ^{x), s,t = 1,..., kn. 


Proposition 2.2. Suppose that the sequence Za satisfies condition D{un{j),kn,ln), where {uniT)}n>i 
is a sequence of real numbers satisfying (1.1) and Bn is a family of subsets of An satisfying the conditions 
of Proposition \K1\ Then, there exists the spatial extremal index, 9 a, if and only if there exists 


lim 

n^+oo 


E 


P 


I \J Z{x) > Unix) 


5 















































and, in this case, we have 


Oa 


lim 


1 


n^+oo f{n)F{Un{T)) 


E ^ 


I \/ Z{x) > Unir) 


Next, we prove that the expected number of exceedances of the level rt„(r) on the blocks Bn’^\ 
s,t = 1,... ,kn with at least one exceedance, converges to the reciprocal of the extremal index 9a- We 
can verify that the greater the clustering tendency of high threshold exceedances (several exceedances on 
the smaller 6a will be. For isolated exceedances of Unin), we have 6a = i- 


Proposition 2.3. Suppose that the sequence satisfies condition D{un{T),kn,ln), where {un('7‘)}n>i 
is a sequence of real numbers satisfying (1-1) and Bn is a family of subsets of An satisfying the conditions 


of Proposition 2.1. If the sequence Za has spatial extremal index, 6a, then 


9A = lim -y 

n-)'+oo 


? E f E I 


{Z(x)>«„(t)} 


^{Z{x)>Un{T)} 

(s,t) 


> 0 


x£B 


If 


lim E I 11 


n^+oo 


{Z(x)>u„{t)} 


\xGB- 


(s,t) 


II{Z{x)>Un{T)} > ^ 1 ^5 


x£B. 


(s,£) 


uniformly in s,t £ {1 ,..., kn}, then 6 a = 1- 


3 Local spatial dependence 


The asymptotic behavior of the maximum of non-stationary and anisotropic random fields, defined 
over discrete subsets of M^, subject to restrictions on the local path behavior of high values is now 
analyzed. 

Criteria are given for the existence and value of the spatial extremal index, which plays a key role 
in determining the cluster sizes and quantifying the strength of dependence between exceedances of high 
levels. To attain this goal, we first introduce a condition for modeling local mild oscillations of the random 
field. This condition is an extension to random fields of the D''(ri„)-condition found in Leadbetter and 
Nandagopalan f[T0] l. 

Throughout this section Bn will denote a family of subsets of An in the conditions of Proposition 2.1. 


Definition 3.1. If V{x) is a finite set of neighbors of a point x £ A and V = {Vix) : x £ A}, then the 
sequence Za verifies condition D”{un, Bn,V) if, as n ^ +oo, 


/c^ sup ^ P 
bIi ' 


Z{x) > Un> y Z{y), \J Z{y) > Un 

y y&V{x) 


0 . 


Although the choice of the family V of neighborhoods can be conditioned by the nature of the practical 
problems under study, here we will illustrate the modeling with a natural choice based on the cardinal 
directions. Therefore, in what follows the initials N, E, S, W will represent, respectively, the cardinal 


6 







directions North, East, South and West. The family of neighborhoods of x € A along directions E and 
N, will be denoted by p,q,r G Z A q < 1, and defined as 

Vg’J = {^(^) : ^ e 

where 


V{x) = {y : (ai(7ri(x)) < 7ri{y) < ap(7ri(x)) A a_g(7r2(x)) < 7r2(y) < a,.(7r2(x))) 
V (7ri(x) = 7ri(y) A ai{TT 2 {x)) < TT 2 {y) < ar{TT 2 {x))) ,x G A}, 

and, for each z G TTi(An), 

... ,a-2i7ri{z)),a-i{TTi{z)),ao{7ri{z)) = 7ri{z),ai{7ri{z)),a2{'n-iiz)),... 

are the points before and after iTi{z), in ascending order. 

In Figure 2 we find an illustration of a neighborhood V(x) G V^'^, p,q,r G Z A q < 1. 

Ihix 


V(x 


ni(^) 


Figure 2: Example of a neighborhood V(x) G with p = 10, q = —1 and r = 11 


The following result proves that, asymptotically, the disjoint events that add to the probability of some 
exceedance of Un over are those where there occurs one exceedance of Un on the location x G Bn’^^ 

and the maximum on the V{x) neighborhood is below Un- That is, regarding these neighborhoods, Z{x) 
is a local maximum. 


Lemma 3.1. Let {un}n>i be a sequence of real numbers and suppose that sequence Za satisfies condition 
D"{un,Bn,V). Then, for each Bn’^'^ G Bn, we have 


P 


I \J Z(x} > Un 


= E ^ 


^Z(x) > Un> 


V ^(y) 

y&V{x) 



As a consequence of Lemma 3.1 the extremal index of Za, 9a, can be viewed, as n —)• + 00 , as the 
mean of tail dependence coefficients of the form 


=-P I V ^iv) > 

y£V{x) 


Z(x) > Unix) , 


(3.4) 
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which are the tail dependence coefficient of Li ([E]). Observe also that if = 1, then X^\v{x),x) 

is the traditional upper tail dependence coefficient introduced far back in the sixties (Sibuya ([E]), Tiago 
de Oliveira ([IS])). 


Proposition 3.1. Let {un('r)}n>i be a sequence of real numbers satisfying (1.1). If sequence Za verifies 
conditions D”{un{T),Bn,V) and D{un{T), kn,ln) then the spatial extremal index of Za, Oa, exists if and 
only if there exists 


Xa 


lim — 

n->-+co j(n) 


Y, AW(P(x),x), 

xeAn 


and, in this case, we have 


9a = 1- Aa- 


Note that some models can verify condition D” only for certain types of neighborhoods. Nevertheless, 
there exist models, as we shall see further on, that verify condition D"{un, Bn, V), for all V. A particular 
case of such models are those that verify a local dependence restriction that leads to isolated exceedances, 
which we shall denominate condition D'{un,Bn) and define as follows: 

Definition 3.2. The sequence Za verifies condition D'{un,Bn) if, as n ^ oo, 

E E P{Z{x) > Un, Z{y) > Un) -^ 0. 


This dependence condition, which bounds the probability of more than one exceedance of Un over 
a block Bn’^^ with approximately elements, together with condition D{un,kn,ln) lead to an unit 


extremal index. In fact, from Proposition 2.1 and condition D'{un{T),Bn) we have 


lim P [ \/ Z{x)<UniT)\ 

n—¥+oo \ '' / 

\a;GAn / 

exp j — lim W p| V 


n^+oo 




exp ( - lim Y Y P {Z{x) > Unix))) + o{l) 

1 n—>-+00 ^^ ^^ 




= exp - lim f{n)P{Z{x) > Unix)) = exp(-r), 

n—>-+00 


which proves the following result. 

Proposition 3.2. Let {un{x)}n'>i be a sequence of real numbers satisfying (1.1). Suppose Za satisfies 
conditions D'{un{x), Bn) and D{un{x),kn,ln)■ Then, 


P I V ^ «n(r) j 

\x£A„ / 


n^+oo 


exp(—r). 


We now present a class of Gaussian random fields that verifies the conditions established in Proposition 
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Example 3.1. Let Z 5 = {Z(x) : x € S} he a standard Gaussian random field on S' C with correlations 
rx,y, x,y £ S, such that 

6 = sup \rx,y\ < 1 , 
x,y&S 

where S = Un>i {^n}n>i is an increasing sequence of sets of isolated points of satisfying 


f{n) 


n^+oo 




fijn) 

'fin) 


n^+oo 


fijn) 

fin) 


>+00 


0 , i = l, 2 , 


with {kn}n>i and {ln}n>i sequences of integer numbers verifying kn +00 and In —>■ + 00 . We will show 
that under the following correlation condition, 


sup \rx,y\ <rn = o 
x,y&An 


1 


log/(n) 


the sequence = {Z{x) : x £ An}n>i verifies conditions D{un{T), kn,ln) and D'{uniT),Bn) with 
Unix) = + bf(n), where an = ( 21 ogn)^/^ and bn = ( 21 ogn)^/^ — i( 21 ogn)“^/^(loglogn + loglTr). 

By Corollary 4.2.9 of Leadbetter et al. (El), we have 


^ P{Z{x) > Unix), Z{y) > Unir)) 





< 

E 

\P 


x,y&B)i"’^^ 


< 

E 

K 







—vlir) \ f(n) , , ,, 

. ^1 I + 

1 + \rx,y\ / ki 


where 4> denotes the standard Gaussian distribution function and K \s a constant depending on 6. 

1 


Now, since sup \rx,y\ <rn = o 
x,y&An 


log/(n) 


^ J]] iL|r^,y|expf- 


, by Lemma 4.3.2 of Leadbetter et al. (i), we obtain 


-nlir) \ , ,2 fin) ^2 




innir)) 

~r ' X,V \ / 


< ^ iL|ra;,j;|exp 

x,y&A„ 

= o(l)> 


-uHt) 


+ L 


^,y\ 


+ oil) 


proving that Z^ verifies condition D'{uniT),Bn) 

Condition D^Unir), kn,ln) follows from Corollary 4.2.4 of Leadbetter et al. (0). 

For other related results concerning Gaussian random fields we refer the readers to Piterbarg (|15j). 
Adler ([I]), Berman ([3]), Choi ([!]) and Pereira f[T4]i. 

Proposition |3.1[ which states that the spatial extremal index 9a is asymptotically equal to the mean 
of tail dependence coefficients is illustrated in the following example with a 1 -dependent random field. 
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Example 3.2. Let {Y(x) : x G Z^} be an independent and identically distributed random field with 
common distribution function Fyly) = exp , y > 0, and define 


Z{x) 


' 9Y{x) 

' V 

, y&E{x) 


7ri(x) = 7r2(x) 
7ri(x) / 7r2(x), 


(3.5) 


where 

E{x) = {y : a_i(7ri(x)) < 7ri(y) < ai(7ri(x)) A a_i(7r2(x)) < 7r2(y) < ai(7r2(x))} . 




200 


150 


100 


50 


0 


Figure 3; Simulation of the random field Z{x) defined in 3.5 


We will calculate the extremal index of sequence = {Z{x) : x G An}n>i, where 

An = {x ^ : —n < 7ri(x) < n, (—n) V (vri(x) — n) < 7r2(x) < (7ri(x) + n) A n} 


and 


f{n) = 3n^ + 3n + 1, fi{n) = 2n + 1, f = 1, 2, 


fijn) _^ 

f (u) n—>-+00 


0 , \jAn=Z‘^, 


n>l 


by two different methods. 

Note that the random field Z = {Z{x) : x G Un>i anisotropic and non-stationary with common 

unit Frchet distribution, F{y) = exp(—y > 0. 

Let Z = {Z{x) : x G Z^} be the associated random field of Z, id est, Z{x), x G Z^, are independent 
and identically distributed random variables having unit Frchet distribution. 

For a sequence of real numbers {u„(T)}n>i verifying (1.1), that is ttn(r) = we have 


lim P \J Z{x) < Un{T) = lim exp 


n^+oo 


KX&An 


fin] 


fin) 


= exp(—r). 


On the other hand, for each Dn = {x G An : 7ri(x) = 7r2(x)} and Unir) = it holds 


P f Y Z{x) > UnirU < fiDnF{UniT)) 

\xGDn / 


277 -U 1 _ 

-^/(n)P(u.(r)) ^ 0, 

jin ) n—>-+00 
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and consequently 


P 



< Unir) 


p 


\J Z{x) < Unir) 


—Dn 


+ o{l) 


p| \J y(x) < Un(r) j +o(l) 
(exp(-u„(r)-V9))^("+^)'+^(”+^^+^ ^ exp 

n^+oo 



so it holds 6a = ^■ 

Let us now consider Proposition |3.1| for the computation of the extremal index of Z^. 

Sequence verifies condition D{uniT),kn,ln), for all sequences of integer numbers {ln}n>i and 
{A:n}n>i satisfying (2.2), since an{ln,Un{T)) = 0 for In > 2. 

Considering V = ^ the family of neighborhoods 

V{x) = {y : (a_i(7r2(x)) < TT 2 {y) < vr 2 (x) A a_i(7ri(x)) < 7ri(y) < 7ri(x)) V 

(7r2(y) = vr 2 (rE) A a_i(7ri(x)) < 7ri(y) < 7ri(x))}, 


sequence Za also verifies (rin(r),V)-condition since, for each x ^ Dn, 


< 


( 


P 


Z(x) > Unix) > y Z{y) 
yev{x) 


P{Y{a) > Un{T),Y{h) > Unix)) 


, \J Z{y)>Un{x) 
yesi^’^Yvix) 

= Fy(M„(r)), 


for pairs of different locations a and b, and consequently 

( 

sup ^ P Z{x) > Unix) > y Ziy), y Z{y) > Unix) 

Bny y&V{x) 

— ^n^'X^^Yi'^nix)) + oil) 

= /(b)( 1-exp(-r/9/(n)))2 = o(l). 


We can then apply Proposition |3.1| to obtain 
6a = 1 - Aa 


1 - lim P ( y Ziy) > Unix) |Z(x) > Unix) 

n—^rxri \ * 

\y&v{x) 

P(y(x) > Unix)) 


E 


1 — lim , / / , / 

n^oo P{Z{x) > Unix)) 

xeB(x)-{(a_i(7ri(x)),a_i(7r2(x)))} ^ ^ ^ 


1 

9' 
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4 Application to max-stable random fields 


Max-stable random fields are very useful models for spatial extremes since under suitable conditions 
they are asymptotically models for maxima of independent replications of random fields. Furthermore, 
all finite dimensional distributions of a max-stable process are multivariate extreme value distributions. 

Within these random fields, it is important to identify dependence among extremes. In particular 
detecting asymptotic independence is fundamental and recently some authors have proposed measures 
of extreme dependence/independence with associated tests. With this in mind we compute, in this 
section, the extremal index of the class of max-stable random fields, as a function of well known extremal 
dependence coefficients found in literature, which will provide immediate estimators for Oa- 

One convenient way to summarize the dependence structure of a max-stable random field Zg' = 
{Z{x) : X £ S}, 5 C with marginal distribution F, is through the extremal coefficient, e/, of 
Schlather and Tawn ([l6]), satisfying 

P [ V Z{x) < y J = y G M, / C 

\xei / 

which measures the extremal dependence between the variables indexed by the set /. Its simple interpre¬ 
tation as the effective number of independent variables indexed in I from which the maximum is drawn 
has led to its use as a dependence measure in a range of practical applications. Another way to access 
the amount of extremal dependence of a random field is through a particular case of the tail dependence 
function introduced in Ferreira and Ferreira ([6]), defined as 


A 


ihlh) 

u 


{yi,y2) 


lim P 

t —^-j-CXD 


V 1 


m 

t 


V > 1 

xGl2 



{yi,y2) e K+, 


(4.6) 


provided the limit exists. 

Note that for max-stable random fields, the limit given in (4.6) always exists. 

The function is a measure of the probability of occurring extreme values for the 

maximum of the variables indexed in a region Ii C given that the maximum of the variables in¬ 
dexed in another region I 2 , with Iif] I 2 = 0, assumes an extreme value too. 

At the unit point, we have 


A 


ihlh) 

u 




xeii 


V 1 

xGh 


which is related with the extremal coefficients of Schlatter and Tawn ([E]), e, in the following way 


A 


(hlh) 

u 


(1) 1) — e/i + e /2 — e/^u/2- 


The next result provides a connection between the dependence structure of the sequence of max- 
stable random fields Z _4 = {Z{x) : x G An}n>i and the limit of the sequence {An^^(I4(x), x)}n>i, with 
^n\v{x),x) defined in (3.4). 

We assume, without loss of generality that for each x G An, n > 1, Z{x) has a unit Frchet distribution. 

Proposition 4.1. Let Za = {Z{x) : x G An}n>i be « sequence of max-stable random fields with unit 
Frchet margins. Then 

lim (1 - Xlf'>{V{x),x)) = ev(x)[J{x} - Wfx)- 
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By combining the previous result with Proposition 3.1 we are able to conclude that if the extremal 
index of exists, then it can be computed from ey(^x){j{x} ~ ^v(x)- We shall denote ey(a;)u{ 3 ;} ~ ^v(x) 
simply by x G An, and name them local extremal indices. 

If we consider the mean of local extremal indices for points on a region of An we obtain a regional 
exremal index, formally defined as follows. 


Definition 4.1. Let Za he a sequence of max-stable random fields with unit Frchet margins. The extremal 
index of Za over a region R C An, with '^R < +oo, is defined as 




fix}- 


xGR 


If Za verifies conditions D{un{T),kn,ln) and D"{un{r),Bn) then, for large n, its spatial extremal 
index, Oa, can be viewed as the mean of local extremal indices, as stated in the following result. 

Proposition 4.2. Suppose that the sequence of max-stable random fields Za has marginal unit Frchet 
distribution and verifies conditions D{un{T),kn,ln) o,nd D''{un{r),Bn)- If 

1 — \!ff\V{x),x) —)• 0|a;|, uniformly in x, 

n^+oo 


then 


Oa = lim V 

n^+oo fin) 

^ ^ xeAn 


We consider in the following example an anisotropic and non-stationary random field to illustrate 


Proposition 4.2 


Example 4.1. Let {Y{x) : x G be an independent and identically distributed random held with 
common distribution function Fyiy) = exp ^ , y > 0 and dehne the random held {Z{x),x G 

(E U F) X Z}, with E = {Ak ; A: G Z} and F = {4k + 3 ; A: G Z}, in the following way 


Z{x) = < 


V Y{y) , xGFxZ 

yeUi(x) 

y Y{y) , X G F X Z, 
^ yeU2{x) 


(4.7) 


where 


and 


Ui{x) = {x, (ai(7ri(x)),ai(7r2(x))), (ai(7ri(x)), a_i(7r2(x)))} 


U 2 {x) = {(a_i(7ri(x)),a_i(7r2(x))), (ai(7ri(x)), 7r2(x)), (a_i(7ri(x)), ai(7r2(x)))}. 


The random held Z = {Z{x) : x G (F U F) x Z} is anisotropic and non-stationary with common 
marginal distribution F(x) = exp(—x“^), x G (F U F) x Z. 

Let us consider 

An = {En U En) X Gn = ({4A: : —n < A < n} U {4fc -|- 3 : —n <k< n}) x {—4n,..., 4n + 3}, 
and consequently, 

fi{n) 


/(n) = 8(2n +1)^, /i(n) = 2(2n + 1), / 2 (n) = 4(2n + 1), 


f (n) n—>-|-oo 


0, An — (F U F) X Z. 


n>l 
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The sequence = {Z{x) : x G An}n>i verifies condition D{un,kn,ln), for all sequences of integer 
numbers {Zn}n>i and {kn}n>i satisfying ( 2 . 2 ), since an{ln,Un) = 0 with In > 2 , as well as condition 
D”(unir), Bn,V) where {un{T)}n>i is a sequence of real numbers satisfying (1.1) and V is the family of 
neighborhoods F(x) = {( 7 ri(x), ai( 7 r 2 (x))), ( 7 ri(rc), a 2 ( 7 r 2 (x)))}. In fact, for each x G {EnXGn)A{FnXGn), 
we have 


P Z{x) > Unix) > V zfe). V Z{y) > Unix) 

\ y&Vix) y^Bi^’*Gv{x) 


P z{x) > Unix) > V Z(!/). V Zfe)>u„(r: 

\ ye{Bi°’'^Gv{x))ATT2{y)>TT2{x) 


+ P Z{x) > Unix) > V Z(s). V Ziu) > Unix] 

\ y&y{p y^iBi"’'^Gv{x))/\TT2{y)=TT2ix) 


< PliUnix))^!^ + Fyinnix))^. 


Therefore, for any family Bn in the conditions of the previously established results, we obtain 


kn sup ^ P Zix) > Unix) > \/ Z(y), y Ziy) > Unix) 


x&Bi^’^'> \ y&V{x) y(zB)^’*'Gv{x) 

^ ifin)Fziunix)))‘^ (/(n)Fy(u„(r)))2 

- - ki -+-e-= 


From the dehnition of spatial extremal index of we have 6a = 5 , since, for Unix) = /(n)/r 
8 ( 2 n + 1 )^/t, we obtain 

P I Y Zix)<Unix)\ = -f’ [ V Zix) <Unix) \ P ( \J Zix)<Unix)\ 


= iF^iunix))) 2 ^ (exp(-r))2. 

n^+oo 
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On the other hand, from Proposition [4]^ 

_ 1 P(Z(x) > M„(r),Z(7ri(a;),ai(7r2(x))) > u„(r) V Z(7ri(x),a2(7r2(x))) >'Un(r)) 

^ /(„) p^z{x) > unir)) 

If X G Eji ^ G<fi: 

P{Z{x) > Un(T),Z{7Ti{x),ai{7r2{x))) > Unjr) V Z(7ri(x),a2(7r2(x))) > Unir)) 

PiZ{x) > Unir)) 

_ PjZjx) > rtn(T),Z(7ri(x),ai(7r2(x))) < ttn(r), Z(7ri(x), a2(7r2(x))) > Un('r)) 

P{Z{x) > Unir)) 

_ P(P(ai(7ri(x)),ai(7r2(x)) > Unir)) 

P{Z{x) > Unir)) 


otherwise, if x G x Gn, 


P{Z{x) > Unir), Z(7ri(x), ai(7r2(x))) > Unix) V Z(7ri(x), a 2 {T^ 2 {x))) > Unir)) 

P{z{x) > Unix)) 

_ PjZjx) > UniT),Z{TTi{x),ai{7r2{x))) > nn(r),Z(7ri(x),a2(7r2(x))) > Unix)) 

PjZjx) > Unjx)) 

PjZjx) > Unjx),ZjTTljx),aijTT 2 jx))) > Unjx),ZjTTljx),a 2 jx 2 jx))) < Unjx)) 

PjZjx) > Unjx)) 

_ PjYja-ijxijx)),aijx 2 jx)) > Unjx)) P(y(a-i(7ri(x)),7r2(x)) > Unjx)) 
PjZjx) > Unjx)) PjZjx) > Unjx)) 

Thus, 


Oa 


1 - lim ( + 1) X 4(2ra + 1) 
n^oo \FzjUnjx)) fjn) 

^ Fyjunjx)) (2n + 1) x 4(2n + 1) ^ 
Fzjunjx)) fjn) ^ 


1 — lim 


fjn)FYjunjx)) , 1 , 2fjn)FYjUnjx)) , 1 

X — X 


1 - 


\fjn)Fzjunjx)) 

h?xi + ^xT = l 

r 2 T 2/2’ 


fjn)Fzjunjx)) 


as expected. 

Although, in practical applications the conditions Djunjx),kn,ln) and D”junjx),Bn) are not easy to 
verify, the results of this section highlight the importance of ia the study of locally occurring large 
observations in clusters. In a region R C An, the smaller the values of 0{a;}) x € R, the greater the 
propensity for clustering. 

Note that, beyond the interpretation of the inverse proportionality between the value of 9a and the 
propensity for clustering, small values of 9^x} indicate a strong dependence between Zjx) and {Zjy),y G 
Vjx)}. 

4.1 Estimation of the spatial extremal index 

The spatial extremal index 9a can be computed from the local extremal indices x G An, and 
the latter are simply the diffrence between extremal coefficients at Vjx) U {x} and Vjx), as previously 
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proved. Several estimators have already been studied in the literature (Krajina (|8]), Beirlant et al. ([2j), 
Schlather and Tawn ([E]), among others) for extremal coefficients. Ferreira and Ferreira ([6]) proposed 
a non-parametric estimator based on the following relation 


E{M{C)) 

1-E{M{C))' 


where M{C) = \J F{Z{x)). 

x&C 


It considers sample means and is defined as 


M{C) 


where M{C) denotes the sample mean, 

k 

i=l y&C 


and F, is the (modihed) empirical distribution function of F, 

1 ^ 

2 = 1 


where Z^'^\y), i = 1,..., k, are independent replications of Z{y). 

With this estimator of the extremal coefficient we propose the following estimator for the local ex¬ 
tremal indices, 

^{x} ^{x}UV{x) ^V{x) W ^ An, (^•8) 


which are consistent, given the consistency of the estimators ^{x}uv{x) ey(3,), proved in Ferreira and 

Ferreira (0). 

From Proposition 4.2 we know that, for large n, 


Oa 


1 

fin) 


XGAn. 


(4.9) 


Therefore, if in (4.9) we replace with its estimator 0|3,| we obtain an estimator for the spatial extremal 
index 6a- 

The finite sample behaviour of the estimator 

= TTT E «{»> 

is analyzed on simulated data from the anisotropic and non-stationary random field {Z(x),x G (FUF) x 
Z} considered in Example 4.1. We simulated 10 times k = 100, 500 and 1000 independent random fields. 

Table 1 shows the mean and mean square error (MSE) of the estimates for k = 100, 500,1000 and 
/(ra), n = 1,10,20,30. 

As we can see from the values reported in Table 1, the estimator 6a has quite a good performance, 
with biases around 0.02 for small values of k and around 0.01 for bigger values of k. The values of n 
considered have a small effect on the bias, nevertheless the variance decreases with n and with k. 
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fin) 



k 




100 

500 

1000 

El.] 

MSE[.] 

El.] 

MSEl.] 

E].] 

MSEl.] 

f(l)=72 

0.5253 

7.6e-4 

0.5140 

2.2e-4 

0.5099 

l.le-4 

f(10)=3528 

0.5243 

5.9e-4 

0.5130 

1.7e-4 

0.5098 

9.5e-5 

f(20)=13448 

0.5237 

5.6e-4 

0.5130 

1.7e-4 

0.5095 

9.0e-5 

f (30)=29768 

0.5237 

5.6e-4 

0.5130 

1.7e-4 

0.5095 

9.1e-5 


Table 1: Estimated mean values and mean square errors of estimator 9a_, for the random field of Example 

4.1. 

5 Conclusion 

In this paper we establish existence criteria for the extremal index of a nonstationary and anisotropic 
random field, defined on M^. Under restrictions on the local path behavior of exceedances, that allow 
clustering of high values, we obtain the extremal index as the limit of a sequence of upper tail dependence 
coefficients. Eor the particular case of max-stable random fields, we prove that the extremal index can 
be obtained as a function of extremal dependence coefficients. Based on this relation we give a simple 
estimator of the extremal index and we analyze its performance with an anisotropic and nonstationary 
1-dependent random field. The simulation study results show the good performance of the proposed 
estimator. 


Appendix A: Proofs for Section 2 

Proof of Lemma 12.11 

If all the subsets and s,t = 1,... ,kn, have less than elements, the result is trivial. On the 
other hand, if some has less than In consecutive elements of 7ri(^n), we can eliminate in the 

family : s,t = 1, .../Cnj, 


since 


P 


= P 


( 


/ 

V > 

\xeUt 


\ 

( \ 


\J Z{x) < Un 


< 


In \/ 


(tt^ Ua) n An) — 


fin) 


f{n)F{un) = o\ — 

i^n 
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and 


n p V z{x)<un -n^ V 


Ur 


kji ( 


s,t 


\x&B. 


(s,t) 


Ur 


^ 1-n^ V 

Vxesi"*’*^ 

/ 

< \J Z{x) > Un 

f{n)F{un) = o(l). 


< 


kn^n \J 


tt (vr^ ^(a 


fin) 


With similar arguments, we can eliminate in the family s,t = l,...A:n| the subsets ti G 

{!,..., kn}, that have less than elements. 

To conclude, let us then assume that all the subsets and s,t = 1,2,..., kn, have more than In 
elements. Start by eliminating in each and respectively, the sets and with the highest 
values In values. The resulting sets s = 1,..., kn, belong to S{Tri{An), In), s = 1,..., kn, and 
t = 1,..., kn, belong to 5(7r2(^n), In)- Now consider 

and note that 


< 


+ 


+ 


P 


P 


P 


Zix) < Un - 

P 

V 

Z{x) < Un 

j s,t 




{ 

( 



Zix) < Un — P 


V 

. . _ 

Z{x) < Un 


V 

( 

\J Z{x)<Un\-WP[ V ^i^)^ 


Up I V Z{x)<Un \ - llP\ V 

2o(l) + kla {In, Un) = o(l), 


Ur 


Ur 




since 


/ 


p 


( 


< 


\J z{x) < Un - P \/ Z{x) < Ur 
fiin) j./ \^/ 


V 2knln^J^f{n)F{Un) = o(l), 

! fin) 
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condition D{un,kn,ln) holds for T^a and 


< 

< 

< 

< 


( 

Up V z{x)< 

/ 


Un - 




Ur, 


s,t 




{s,t) 


Ur,, - 




(s,t) 


V zw< 

" * \xeBi^'*^ 

FiUn) 

S t 

lnF{Un) {knfl{n) + knf 2 {n)) 

V 2knln^^^p\f{n)F{Un) = o(l). 

: f{n) 


P I y Z(x) < Ur 


Proof of Lemma 12.21 

(s t) 

First, we prove that the disjoint subsets Bn ’ ^ s,t = 1 , ...,kn, constructed through the method mentioned 
before Lemma 2.2, has approximately elements. 

1. In 'Ki{An), let us consider the set of the first elements maximally constructed such that 

P{Z{x) > Un) < Y P{Z{F > '^n)- 

In 7ri(74„) — iF let us consider the maximal set, of the first elements such that 

P{Z{x) > Un) <Y P{Z{x) > Un). 

a;G7r^^(7f^)) P An 

Similarly, we obtain kn disjoint subsets of An^ 

TT^iFF) r^An,.. . n An. 

2. For each one of the previous subsets, let us consider an analogous decomposition using projection 
T^2. 

In Tr 2 (An), we consider the set that consists of the first elements maximally chosen such that 

P(Z(x) > Un) < ^ y] P(Z(x) > Un). 
xg7r^^(j(®a)) p|7r“l(/(s)) p ^ x£A„ 

Using the same technique, we obtain the following kn subsets of 7rj“^(/(^)) n An, 

BF^) = n n An, t = l,...,kn. 


Next, we will prove that ~ Since 


Y P{Z{x) > Un) < 

x£tt~^{J(‘‘’*)) n 7 rl“^(/('>)) n An 


L2 Y1 Pi^iF > Un) 
” X&An 
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we obtain 


< 4/(n), 

and, from the maximality criteria used in the construction of it follows that 


+ o ( > ^f{n)F{un). 


Therefore, 


kn^n J 


kl ( 1 

1 > . . > 1 - 


^ - F(n) \knln)' 

/C„ 


which allows us to conclude that ~ since nF^oo 


Now, 


P I Y Z{x) <Unj 

\x£An ) 


- P 


( 

\J Z{x) < Ur. 


< P{Z{x) > 

xeAn- Us,tB^’^'^ 


Ur. 


An 1^ Pn 

V S,t 


{s,t) 


-f{n)F{ur. 


and 


f{n) 




s,t 


s,t 


fin) 

1 v-tJP: 


{s,t) 


1 — — < 1 - A 

kl ^ /w - 


/(’^) 


ttp. 


(s,t) 


fH 

“fcH“ 


= o(l), 


which concludes the proof. 


Proof of Proposition |2.2| 

From Proposition |2.H we have 


lim P [ \/ Z{x) < Unix) ] 

\xGAr. J 

lim n ^ I V Z{x)<Unii 


n^+oo 




n^+oo 


n^+oo 


tion of 9 a- 


E 

- In 1 - P 


\ V 

/ 

E 

p[ y z{ 


VxeBl'’’*' 

’eB„ ^ ( 



.XdB. 


(s,t) 


20 















Proof of Proposition |2.3| 

Follows from Proposition |2.2[ since 


6 a = lim 


X] V Z{x)>Un{T) 




n^+oo f{n)F{Un{T)) 

P j \J Z{x) > Unix) 


= lim -y 

n->-+oo 


E 


\x£B, 


(s,t) 




^F{Un{T)) 


P I ^ ^{Z{x)>u„{t)} > 0 


= li™ 7*2 X] ' 

n-s>+oo 


\x^B. 


(s,t) 


® E “ 


{Z(x)>Un{T)} 


\xGB, 


(s,t) 


Appendix B: Proofs for Section 3 

Proof of Lemma 13.11 


We can write 


p V 


'^n 


\xGB. 


(s,t) 


/ 


= p 


u 


\xeBl 


s,t) 


E ^ 




xGB. 


(s,t) 


Z{x) > Un> y Z{y) 

:7ri (a:)>7ri (j/)V(7ri {y)=-n;i (x)Att 2 {y)>TT2 (a;))} 

Z{x) > Un> y Z{y) 

{y&Bi^’*^:TTi{x)>wi{y)W{TTi{y)=wi{x)ATr2(y)>TT2{x))} 


^ piz{x)>Un> y Z(y)|+I 


x£B. 


(s,t) 


y&V{x) 


since, by D"{un-,Bn,V^E% )-condition, we have 


x£B 


P Z(x) >Un> y z{y) 

(“B y y&V{x) 

I 


E P 


x£B- 


(s,t) 


Z{x) > Un> 


\ 


V 


ziy) 


:7ri (x)>7ri (j/) V(7ri {y)=Tri (x)Att 2 (y)>TT2 (a;))} 


/ 


< 




sup Y P Z{x) >Un> y z{y), y Z{y) 


\ 


y&v(x) y(ZB^^A_V{x) 


> Un\ = O I Vs 
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Remark 1. Instead of the decomposition 0/{\/^g^(s.t) Z{x) > Un} as the union of events {Z{x) > Un}, 

(s t) ^ 

X G Bn’ , where x is the location with the highest coordinates, previously considered, we can consider 
other decompositions. For example, the events {Z{x) > Un} where x G is the location with the 

lowest abscissa and biggest ordinate, lead to the following decomposition 


V 


> Un} = 


U 


) x£B, 

Now, denoting by p,q,r & Z, the family of neighborhoods 

V{x) = {y : {a-p{'Ki{x)) < -Ki{y) < ai(7ri(x)) A a_q(7r2(x)) < 7r2(y) < ar{'K 2 {x))) 

V(vri(y) = 7ri(x) A ai{'K 2 {y)) < T^ 2 {y) < ar{TT 2 {y))) 

we obtain an analogous result to Lemma 3.1, if we assume condition D”{un, Bn,V^^). 

In fact, we can decompose the event W!Z{x) > Unir)} in eight different ways corresponding 
to the different forms of starting from x and getting around T?, along the directions {{ak{7ri{x)), 0), k G 
{-1,1}} and {(0,afc(7r2(x))), A: G {-1,1}}. 




Z{x) > Un> y z{y) 

: 7 ri (y) <tti (x)V(Tri (y)=ni {x)/\-K 2 {y)>T^2 (a;))} 


. 


Proof of Proposition |3.1| 

and the definition of the tail dependence coefficient \y{V{x),x), it follows that 


From Lemma 


3.1 


\l Z{x) > Unix) 


\x£B. 


(s,t) 


P I Z{x) > Unix) > y Z{y)\+o(^^ 

FiUnix)) X] ^ f V Z{y) <Un{x) \ z{x) > Un{x)\ +o(^'] 
\v&V{x) j ^ 

F{Un) Y (1 “ +0 fiV 
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Thus, 


lim P [ \/ Z{x)<Un{T)] 

/ 

= exp ( - lim ^ pi y Z{x) > Unir) 

I n—¥oo \ 

lim (p(un(r)) Y (l-^n^^(^(3;)>3;)) + or 

n^oo I \ 

= exp I - lim P(M„(r)) Y ^ (1 “ 

' n—>-00 z_—/ 


= exp 




= exp - lim F{un{j)) V (1 - aM(F(x),x)) 

\ n—>-00 ^ 

V x&An / 


= exp —T 


which concludes the proof. 


) 


Appendix C: Proofs for Section 4 

Proof of Proposition |4.1| 

Since Unir) = and F{y) = exp(—y“^), we have 

lim aP(P(x),x) 


lim P ( V Z{y) > ® 


\y&v{x) 


z{y)> 


fin) 


itapM V z(rf)>p(fP) 


yy&V{x) 


F{Z{x)) > F 


fin) 


lim P [ \/ FiZiy)) > exp f - 

1^00 I ^ \ 


. yGV(a:) 


fin) 


F{Z{x)) > exp ( — 


fin) 


= A 


lim P V F{Z{y)) > 1 - T^T 

\ ^ fi^) 

\y£V{x) 

{V{x)\{x}) 


FiZix)) > 1 - 


fin) 


u 


f T ^y(a:) ^V{x)[J{x} 
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Proof of Proposition |4.2| 

Results from the fact that {1 — Ai^^(P(x), x)}n>i converges to uniformly in x, since this implies 
that for every e > 0, there exists a natural number m = m{e) such that, for every n > m, 


— j;(i-A(:)(p(x),x)) < 

xGA„ 


1 

fin) 


Y1 i^{^} + 

X&An 



+ e 


and 


.ti. 


1 

fin) 


X&An 


e) 



— e. 
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